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ABSTRACT 

We report on a 120 ks XMM-Newton observation of the galaxy cluster Abell 2597. Results 
from both the European Photon Imaging Camera (EPIC) and the Reflection Grating Spec- 
trometer (RGS) are presented. From EPIC we obtain radial profiles of temperature, density 
and abundance, and use these to derive cooling time and entropy. We illustrate corrections to 
these profiles for projection and point spread function (PSF) effects. At the spatial resolution 
available to XMM-Newton, the temperature declines by around a factor of two in the central 
ISOkpc or so in radius, and the abundance increases from about one-fifth to over one-half 
solar The cooling time is less than lOGyr inside a radius of 130kpc. EPIC fits to the central 
region are consistent with a cooling flow of around 100 solar masses per year Broad-band fits 
to the RGS spectra extracted from the central 2 arcmin are also consistent with a cooling flow 
of the same magnitude; with a preferred low-temperature cut-off of essentially zero. The data 
appear to suggest (albeit at low significance levels below formal detection limits) the presence 
of the important thermometer lines from Fe XVII at 15, 17 A rest wavelength, characteristic of 
gas at temperatures ^ 0.3 keV. The measured flux in each line is converted to a mass deposi- 
tion estimate by comparison with a classical cooling flow model, and once again values at the 
level of 100 solar masses per year are obtained. These mass deposition rates, whilst lower than 
those of previous generations of X-ray observatories, are consistent with those obtained from 
UV data for this object. This raises the possibility of a classical cooling flow, at the level of 
around 100 solar masses per year, cooling from 4keV by more than two orders of magnitude 
in temperature. 
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1 INTRODUCTION 

Abell 2597 is a relatively nearby (z ~ 0.08), Abell richness class 
cluster of galaxies. Rich in features, it has been studied extensively 
in many wavebands. 

The central cD gala xy contains the radio sourc e PMN J2323- 
1207 (PKS B2322-12) jWright & Otrupcekl 1199(1 iGriffith et all 
11994. The total flux density of this sourc e is ~ 1.7Jy at 1.5 GHz 
/McNamara et al.'200lV Owen et alll99l . and ~ 0.3 Jy at 5.0GHz 
(Ball et al. 1993), implying a fairly steep spectral index a = — 1.3 
{Sy t>= v"). The source is physically small, ~ 7 x 5,5 x 2arcsec 
at 1.5, 5.0GHz respectively. The source is only moderately re- 
solved at these frequencies, but at the higher frequency has an elon- 
gated appearance suggesting a double-lobed structure. At 8.44GHz 
iSara zin et al. 1995), the source is resolved into a nucleus and 
two lobes (~ 5 arcsec across), orientated roughly NE-SW. A jet 
is clearly seen leading from the nucleus on the SW. This jet un- 
dergoes a sharp (~ 90°) bend about 0.5 arcsec from the nucleus. 
Assuming a typical spectral index for the jets of —0.7, the inferred 
lobe spectral index is ~ —1.5, suggestive of confinement and syn- 
chrotron ageing. Very Long Baseline Array (VLBA) 1.3, 5.0GHz 
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observations of the central 0.3 arcsec <TavloretalJll999t) show an 
inverted spectrum (a = 0.6) core, together with straight, symmetric 
jets emanating from both sides. 

Abell 2597 has previously been observed in X-rays by 
variou s observ atories, for examp l e: Einstein iCrawford et'all 
11989^; ROSAT I'Sarazin et al." 1995; 'Sarazin & McNamara"l997l 
fPeres et al. 1998); A5C AJ White 2000); C/ianrfra I McNamara e t al] 
2001); XMM-Newton Istill & Mushotzkvl 120021) . Traditior ially, it 
has been classified as a moderately strong cooling flow (e.g. lpabianl 
1994) cluster, with inferre d X-ray mass depo sition rate ( Mq yr^'): 
270 ±41 (ROSAT PSPC: fPeres erai]ll998l r^^oX = 152+^^ kpc); 
259^178 (A5CA;iWhiia2Q03). 

A short (~ 20ks of good time) Chandra observation 
iMcNamara et al. 2001) highlighted the interaction between the ra- 
dio source and the X-ray gas, through the presence of so-called 
'ghost cavities' of low surface brightness, seemingly coincident 
with spurs of low frequency (1.4GHz) radio emission. Such bright- 
ness depressions are believed to be buoyantly rising bubbles (e.g. 
_Blrzan et al. 2004) associated with a prior (~ 10'^ yr previously, in 
this case) outburst of the central radio source. The small physical 
size of the radio source means, however, that these issues are not 
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amenable to study with XMM-Newton, given its ~ 5 arcsec point 
spread function (PSF). 

In the Chandra and XMM-Newton era of high spatial and 
spectral resolution, the cooling flow model is undergoing re- 
interpretation. The cause is quite simply an absence (or at best 
a severe weakness) of the spectral features expected from gas at 
the low-end (< IkeV, say; though there does not appear to be a 
fixed absolute cut-off) of the X-ray temperature regime in the ob- 
served spectra of cluster c entral regions. Red uctions in temperature 
by factors of a few (e.g. Kaastra et al."2004') with radius are seen. 
The expected emission lines for several important low-temperature 
species (e.g. the FeXVII 15 and 17 A lines) do not seem to be 
present at the levels predicted by simple models. This is a trend 
that appears to be repeated in many clusters that have traditionally 
been thought to harbour cooling flows (e.g. Peterson et al. 2003). 
Many explanations (e.g. .Fabian et aL ,2001. : .Peterson et alj.2001i') 
for this effect have been put forward; some involving suppres- 
sion of the gas cooling, some involving departures from standard 
coUisionally-ionized radiative cool ing. After a flirtation with ther- 
mal conduction (but see e.g. ISoker 2003), theoretical studies at 
present seem to be concentrat ed on some form of distributed heat- 
ing from a central AGN (e.g. IVecchia et aljboO^I ; [kevnolds et al. 
l2005t iRuszkowski et alj|2004 . effected by buoyant plasma bub- 
bles. 

Observations at ultra-violet (UV) wavelengths suggest that 
Abell 2597 may be a particularly interesting object for study in light 
of these issues. Abell 2597, together with the z = 0.06 cluster Abell 
1795, was observed with the Far Ultra violet Spectros copic Explorer 
(FUSE-. Moos et al. 2000) mission bv lOegerle etall 1.2001). FUSE 
is sensitive to the Ovi 1032, 1038 A resonance lines, which are 
strong diagnostics of thermally radiating gas at temperatures ~ 
3 X 10^ K. In a traditional cooling flow model, in which the ulti- 
mate fate of the hot gas is to cool to very low temperatures, strong 
UV emission is expected in such lines as gas cools out of the X-ray 
temperature regime; thereby potentially forming a bridge between 
the hot. X-ray emitting gas and the cool. Ha (see below) radiat- 
ing gas. O VI 1032 A emission was detected in Abell 2597, with an 
inferre d luminosity 3.6 x l O'^'^ergs^'. Using simple cooling-flow 
theorv. lOegerle et alJ llOOlh calculate the associated mass flow-rate 
as 40 M0 yr^^ (within the effective radius of the FUSE aperture at 
this redshift, ~ 40kpc), for an intermediate case between the limits 
of isobaric and isochoric coolin g. The ROSAT HRI mass deposi- 
tion rate of ISarazin et alj fi995l) . ~ 330 M0yr^' within a cooling 
radius (defined by feool < lOGyr) ~ 130kpc, is around three times 
larger (applying a simple linear M r scaling). The weakness of 
cluster cooling flows in comparison to the predictions of the tra- 
ditional model is thus again demonstrated, but extended down to 
the UV regime. This under-luminosity at UV temperatures may be 
contrasted with the situation at Ho: temperatures, where the Abell 
2597 emission nebula (see below), like other such nebulae, is over- 
luminous, eve n when com pared to classical, high X-ray mass de- 
position rates iVoit & Donahue, 199 7). 

In contrast, Abell 1795 was not detected in Ovi. Similarly, 
iLecavelier des Etanes et all i2004h . also using FUSE, were unable 
to detect O VI emission from either of the low redshift (z ~ 0.07), 
low Galactic column (N^i ~ 4 x 10^" cm^^), rich clusters Abell 
2029 and Abell 3112. These results highlight the apparently un- 
usual (within the currently favoured non-cooling flow paradigm) 
nature of Abell 2597, and motivate a more detailed X-ray study of 
its properties. 

The central (~ 15 arcsec in diameter) regions of Abell 2597 
harbour an optical emission line nebula, as is often (and exclu- 



sively) the case in cooling flow clusters. The Ha + [A'll] luminos- 
ity is 2.7 X 1 0^^ergs~', making it among the most luminous in 
the sample of lHeckman et alJ il989l) . The Ha luminosity is many 
(~ 300) times greater than that expected from extrapolating the 
classical ~ 10^ KX -ray mass flow rates down to ~ lO'* K, a general 
feature of such nebulae. The discrepancy is made even worse if the 
relative sizes of the nebulae and X-ray c ooling regions ar e taken 
into account. The deep optical spectra of IVoit & Donahu3 il997h 
reveal that the Abell 2597 nebula has T ~ lO'^K, Z ~ 0.5 Zq, 
/le ~ 200 cm^-', and is significantly reddened by (intrinsic, owing 
to the low Galactic column) dust. The low inferred column density 
^HII ~ 3 X lO'^cm^^ relative to that of Hi suggests that thin ion- 
ized H II layers surround cold neutral H I cores. The nebula has a 
low ionization parameter, most species being only singly ionized. 
Photoionization from hot stars is the only ionization mechanism 
not ruled out by the observations, but even the hottest O stars (in 
isolation) would not be able to heat the nebula to the observed tem- 
perature. Some form of heat transfer from the intracluster medium 
(ICM) may supply a comparable amount of energy to the nebula. 

The co re of the central cD exhibits a s i gnificant blue optica l 
excess (e.g. iMcNamara & 0'Comielllll993t ISarazin et ai]ll995h . 
HST observations iKoekemoer et alll999h resolve this excess into 
continuum emission around and across the radio lobes, and sev- 
eral knots to the south-west. Line-dominated (thereby not optical 
synchrotron or inverse Compton in origin) excess appears in bright 
arcs around the edge of the radio lobes, and more diffusely across 
the lobe faces. Such geometry is suggestive of an emission shell 
surrounding the radio lobes. A bright optical knot is also seen co- 
incident with the southern radio hot-spot, consistent with a massive 
clump responsible for the deflection of the southern jet. The cen- 
tral region also shows strong dust obscuration, resolved into several 
large (several hundred parsecs) regions. 

In the infra-red, vibrational molecular hydrogen emission 
lines characteristic of ~ 1 000-2000 K collisionally ionized gas 
have been detected from the irmermost ~ 3 arcsec of the central 
cluster galaxy I Falcke et alll998l;lDonahue et all 2o'oo';'Ed ge etall 
^002). As with the optical nebula, there is too much emission as 
compared to the expectation from even an inordinately massive 
classical cooling flow, and yet too little mass for this to be the final 
reservoir of a standard, long-lived flow. Furthermore, the molecular 
gas is dusty, and therefore unlikely to be the direct condensate of 
the hot ICM (where dust is rapidly destroyed). Again as with the 
optical light, there is a complex filamentary structure that seems 
to be spatially associated with the radio source. The optical and 
infra-red properties are both consist ent with UV heating o f the gas 
by a population of young, hot stars (Donahue et al. 2000) (SFR ~ 
few M p) yr^' ). Abell 2597 also has a possible CO detection 

(Ed3 

I2OOII) ; indicative of very cold, molecular gas. 

The Galactic column density in the direction of Abell 
2597 (/ = 65.4, b = -64.8) has the relatively low value 2.5 x 
10^°cm~^ ( Dickev & Lockman'l990; Stark et al. 1992). Recently, 
iBames & Nulsen ^2Q03) have shown that the 5a limit on any fluc- 
tuation in the foreground Hi column in this direction, on scales 
~ 1 arcmin, is 0.43 x 10^" cm^^ (17 per cent). 

Intrinsic H I 21 cm absorp tion towards the central radio source 
was detected by lO'Dea et alj il994 . A narrow (~ 200km s"'), 
redshifted (~ 300kms^'), spatially unresolved component coin- 
cident with the nucleus is inferred to be a gas clump falling onto 
the cD. A broad (~ 400kms^'), spatially extended (~ 3arcsec, 
i.e. the size of the radio source) component at the systemic ve- 
locity is consistent with being associated with the (hence photon 
bounded) Ha. The inferred column densities and masses for the 
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narrow and broad components are A'h ~ {8.2,4.5} x lO^^rsCm^^; 
M ~ {3.5,7} X 10^ Ts Mq; with Ts the spin temperature in units of 
lOOK. Int rinsic, redshifted H I was also seen in the VLB A obser- 
vations of iTavlor et alj Jl999h . and interpreted as a turbulent, in- 
wardly streaming, atomic torus (scale height < 20pc) centred on 
the nucleus. 

Our chosen cosmology has Hq = 50kms^' Mpc^', ilj^ = 

I. 0, and £1\ = 0.0. Adopting a redshift z = 0.0822 (see Sec- 
tion l5.2.ll . the luminosity distance is Dl = 503 Mpc, the angular 
diameter distance is Dp^ = 429 Mpc, and the angular length scale is 
2.08 kpc arcsec^ ' . Data processing was carried out using the XMM- 
Newton Science Analysis System' (SAS) version 5.4.1 (aka xmm- 
sas_20030 110.1802). Spectral analyses employed XSPEC^ version 

II. 3. Unless otherwise stated, quoted uncertainties are IcT (68 per 
cent) on a single parameter. 



2 EPIC DATA REDUCTION 

Abell 2597 was observed by XMM-Newton on 2003 June 27-28 
(observation ID 0147330101), during revolution 0650. The nom- 
inal length of the observation was 120ks, unfortunately the ex- 
posure was heavily affected by background flaring (see Fig. 0. 
The thin optical blocking filter was applied to the European Pho- 
ton Imaging Camera (EPIC). The MOS (Turner et al. 2001) and pn 
iStnider et alJl200lh detectors were both operated in Full Frame 
mode. 

EPIC events files were regenerated from the Observation Data 
Files (ODFs) using the SAS meta-tasks epchain and emchain. 
For emchain, standard options were used, except that bad pixel de- 
tection was switched on. epchain was run in two steps in order to 
merge in the instrument house-keeping Good Time Interval (GTI) 
data (which are otherwise ignored) before the final events list was 
created. 



2.1 Background filtering and subtraction 

The EPIC background is comprised of a number of components: 
proton, particle, photon. Soft protons in the Earth's magnetosphere 
scatter through the mirror system, and are responsible for periods of 
intense background flaring, where the instrument count rate can rise 
by orders of magnitude. Such intervals can only be dealt with by ex- 
cising them from further analysis using a rate-based GTI filtering. 
Energetic charged particles that pass through the detector excite flu- 
orescent emission lines in several of the constituent components. 
The strong spatial variation of some components of this e mission 
jprevberg et id]l2002al Iprevberg et aljE002bl lLumbll2002l) makes 
it preferable to extract a background from the same region of the 
detector as the source. For the analysis of extended sources, the use 
of blank-sky background templates l ,Lumb»2002) is therefore to be 
desired. 

For Abell 2597, cluster emission is detected up to around 8 
and 9keV in the MOS and pn respectively. The pn exhibits strong 
fluorescence lines of Ni, Cu and Zn around 8keV. The only non- 
negligible fluorescent line visible in the pn data below 8keV is the 
1.5keV Al line. Unlike the Cu (fo r example) line, this l ine does 
not vary spatially over the detector I Frevb erg et all2002lj) . since it 
originates in the camera housing. Cutting off the spectra at 7.3 keV 
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Figure 1. MOS2 light curve formed from 100 .s time bins. Upper panel: total 
exposure, with rejected periods of high count-rate in grey, accepted periods 
in black. Lower panel: only the selected Good Time Interval (GTI) is shown. 



removes very few counts, and means that the only background line 
remaining in the pn is not spatially varying. This frees us from the 
need to extract background spectra from the same spatial region of 
the pn detector, i.e. enables us to use a local background. This is 
desirable because it removes the vagaries of scaling and varying 
cosmic background necessarily associated with the use of a blank- 
sky background (see below); and it allows us greater freedom in the 
filtering of background flares. Additionally, comparison of the off- 
source pn spectra for the Abell 2597 and blank-sky fields reveals 
that the 8keV fluorescence lines are shifted to a lower energy by 
about 20 eV in the latter. This almost certainly represents a change 
in calibration between the processing of the blank-sky templates 
and our re-processing of the ODFs for Abell 2597 (the effect of 
such a shift would of course be largely restricted to an increase in 
in those channels lying around the fluorescent line energies). 

The only significant fluorescent lines in the MOS data are the 
Al and Si lines at 1.5 and 1.7keV respectively (the Cr, Mn, and Fe 
lines around 6keV are seen to be negligible). Both of these features 
exhibit strong spatial variation (Lumb 2002) over the MOS field of 
view (hereinafter fov). Using a local background for the MOS is 
therefore not as viable a proposition as it is for the pn. 

In order to be a ble to make u se of the blank-sky EPIC back- 
ground templates of 'Lumb' we use the same form of fil- 
tering process as employed in the construction of those files. That 
is, we make a light-curve (e.g. Fig. in 100 s time bins of the 
high energy (PI > 10000) single pixel (PATTERN == 0) events 
across the whole detector. The standard XMM event selection flags 
#XMMEA_EM and #XMMEA_EP are applied for the MOS and pn respec- 
tively (these exclude cosmic rays, events in bad frames, etc.). 

Our preferred filtering method is to form a histogram of the 
number of events in the 100 s time-bins of the light curve, fit the 
core of the distribution (i.e. ignoring the high count-rate tail) with a 
Gaussian, and exclude periods where the count-rate lies more than 
na away from the mean. Unfortunately, for MOS 1 (and the pn, but 
this is not an issue since we use a local background in this case) the 
2cT count-rate thresholds lie above the filtration limits used in the 
creation of the background events templates (0.2 and 0.45 cts^' for 
the MOS and pn respectively). In order to make use of the MOS 1 
background template we therefore apply the same absolute count- 
rate filtering to this instrument as was used to make the background 
events file. For M0S2, the 2a count-rate threshold is 0. 13cts^', 



4 R. G. Morris and A. C. Fabian 



and this absolute value is used to filter both the observation events 
file, and to re-filter the background template events file to the same 
level. For the pn, the 2a count-rate threshold is 0.54cts^' . 

We strengthen the filtering by requiring that all GTIs be at 
least 10 minutes long. The GTI files are aligned on the frame 
boundaries of the relevant events files before being applied. The 
mean remaining ONTIME of the events files after filtering in this 
way are 39 (pn), 37 (MOSl), and 38 ks (M0S2). 

It is known that the particle background can vary from 
observation to observation at around the ten per cent level 
jprevberg et alij,2002a) . We therefore scale the products (images 
and spectra) of the blank-sky MOS events file before they are 
used to correct the source products. The scaling factor is ob- 
tained from the ratio of the count-rates outside the fov, where 
the only events should be the result of particle-induced fluores- 
cence rather than from actual photons passing through the mir- 
ror system. In detail, we find the number of counts remaining 
in the GTI-filtered data after applying the evselect filtering 
criterion #XMMEA_16 && FLAG k 0x766a0f 63 == && PI in 
[200:12000] && PATTERN <= 12, where the flag expression 
selects events outside the MOS fov but with no other non-zero flag 
settings. Combining this figure with the exposure time gives an out- 
of-fov count-rate which is used to normalize the background tem- 
plate products. In this way, we obtain scaling factors of 1.02 and 
0.91 for the MOSl and MOS2 backgrounds respectively. To al- 
low for some uncertainty in the normalization, we assign a 10 per 
cent systematic error to the blank-sky products, which is added in 
quadrature to the statistical errors. 

The SAS task attcalc was used to reproject the blank-sky 
templates to the same sky position as the Abell 2597 observation, so 
that the same sky coordinate selection expressions could be applied 
to both the science and background fields. 

B elow 5keV o r so, the cosmic X-ray background (CXB) dom- 
inates <Lumbll2002h over the internal camera background. If there 
were no significant variation of the CXB over the sky, then the 
blank-sky templates would also represent this background compo- 
nent. To allow for any changes in the CXB in the region of the 
target, a so-called double background subtraction approach can be 
used. The large field of view of XMM-Newton, together with the 
redshift of Abell 2597, result in significant areas of the detectors 
being essentially free from cluster emission. This enables us to ex- 
amine the background spectrum from such regions and compare 
it with those of the (scaled) blank-sky templates extracted from 
the same detector region. For this purpose we extracted spectra 
from an annulus lying between 7 and lOarcmin (avoiding regions 
near the edge of the fov) from the cluster centre, which we take 
to be located at the emission peak, in sky coordinates {X,Y) = 
(23720.5, 23720.5). Non-ICM sources in this region were excluded 
with circular masks of radius at least 30arcsec. The BACKSCAL 
FITS header keys were corrected for the area so removed. 

Initial comparison of the spectra extracted from this region 
of the Abell 2597 field reveals that the agreement with the blank- 
sky spectra scaled according to the out-of-fov count-rates is good; 
despite the relatively poor quality of the science observation as a 
whole (i.e. the heavy background flaring). There is something of a 
slight excess of soft (< 2keV) counts in the Abell 2597 field, but 
not to any great degree. 

In summary, we use a local (extracted from 7~10arcmin in 
radius) background for the pn. For the MOS, we use the blank-sky 
templates to deal with the spatially varying Si and Al background 
lines, corrected by a double-background subtraction approach to 
account for any (minor) variance in CXB between the science and 



background fields. For the MOS, we make use of single, double, 
triple, and quadruple pixel events (i.e. PATTERN ^ 12). For the 
pn, we use single and double events (PATTERN ^ 4). In both cases 
only events with FLAG == are considered. 

2.2 Out of time events 

When the pn is operated in Full Frame mode, as is the case here, 
bright sources can be significantly affected by so-called out of time 
(OOT) events. These are caused by photons which arrive while the 
CCD is being read out (charge shifted along columns towards the 
readout node). Such events are assigned a wrong RAWY coordi- 
nate and hence an improper CTI (charge transfer inefficiency) cor- 
rection. In an uncorrected image, the OOT events due to a strong 
source are visible as a bright streak smeared out in RAWY. OOT 
events act to broaden spectral features. 

OOT events can only be corrected for in a statistical sense. 
We follow the standard procedure of re-running the SAS meta- 
task epchain with the option withoutof time=y passed to the 
epevents task. This generates an events file that consists entirely 
of OOT events, distributed over the range of RAWY present in the 
data. OOT images and spectra can then be produced in the normal 
way, and scaled and subtracted from source images and spectra to 
correct them for the OOT events. 

The appropriate scale factor is obtained from the expected 
fraction of OOT events, which depends on the ratio of the read- 
out time to the frame time. For the pn in Full F rame mode, this is 
4.6/73. 3ms = 6.3 per cent jStriider et ai]|200lL table 1). A shght 
correction to this factor is needed because the version of epevents 
used does not take account of the fact that in pn Full Frame mode, 
the first 12 (of 200) rows are marked as bad. Accordingly, the cor- 
rect scaling factor is (1 — 12/200) *6.3 = 5.9 per cent. In practice, 
OOT events make little difference to the pn spectral fits. 



2.3 Response files 

The lower limit for the spect ral fits was p laced at 500 eV for the 
MOS (following the advice of lKirschl2003l for post-revolution 450 
data; in any case the MOS results are relatively insensitive to the 
adopted low-energy cut-off). For consistency (and to avoid low- 
energy residuals) the same limit was applied to the pn data. The 
high energy cut-off was placed at 7.3 keV, slightly below the maxi- 
mum detected energy from the cluster, to avoid the background pn 
fluorescence line complex at 8keV (as described in Section lTTl . 

Spectra were grouped to a minimum of 20 counts per channel 
to ensure applicability of the statistic. 

Redistribution matrix files (RMFs) and ancillary response files 
(ARFs) were generated using the SAS tasks rmf gen and arfgen 
respectively. RMFs were made using the standard SAS parameters. 
ARFs were made in extended- source mode, using a low-resolution 
detector-coordinate image of the spectral region. The image pixel 
size was chosen according to the size of the extraction region 
so as to adequately sample the variation in emission, bearing in 
mind that the version of arfgen employed does not model vari- 
ations on scales < larcmin. Strictly sp eaking, the images shoul d 
be exposure-weighted (i.e. de-vignetted; ISaxton & Siddiauil200'3) . 
but we find this makes almost no difference in practice in this 
case. Bad pixel locations were taken from the relevant events file. 
The approach of generating ARFs appropriate to the emission pat- 
tern in the spect ral extraction region differs from the 'weight- 
ing' method (e.g. lAmaud et alj|200ll) more frequently employed 
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Table 1. Global (l^arcmin) spectral properties. 




Figure 2. pn count image of Abell 2597, selected with FLAG == kk 
PATTERN ^ 4 && PI in [200:15000]. The image has been GTI fil- 
tered and conected for OOT events as described in Section |2| but no 
vignetting or background correction has been applied. Image pixels are 
4arcsec across, roughly the size of the PSF core. 

in XMM-Newton spectral analyses. We find that the differences in 
the temperature p rofiles etc. obtained usi ng the two approaches are 
negligible (as per lGastaldello et al^2003^ . 



3 GLOBAL PROPERTIES 

The pn image of the cluster is shown in Fi g. ID As in the ROSAT 
observation of iSarazin feMcNamarljiggTh . there are a number of 
non-ICM sources in the field of view, though only two of these 
are within 6arcmin of the cluster centre. The cluster emission itself 
displays a relaxed, somewhat elliptical form. 

In order to examine the global spectral properties of the clus- 
ter, and to check the agreement between the three EPIC detectors, 
we extracted spectra from an annulus of radius l-4arcmin centred 
on the cluster emission peak (12000: RA = 23h25ml9.8s, Dec = 
-12d07m26s; sky coordinates: X = 23720.5, Y = 23720.5). This ra- 
dial range was selected to excise any cooling flow region and so 
that the background contribution was insignificant, in order to have 
the most simple spectral shape possible. The results of basic single 
temperature fits are presented in TableQ 

When a fixed. Galactic absorption is used (model A), the pn 
has a temperature somewhat below that of the MOS instruments, 
and a higher reduced x'^- If the absorbing column is allowed to vary 
(model B), then all three detectors prefer a value significantly be- 
low Galactic, with the pn Nn being much lower than that of either 
MOS. Importantly, however, freeing A'h results in the MOS and 
pn T and Z values coming in to much closer agreement with each 
other. The conclusion we may draw from this is that all three detec- 
tors exhibit some form of soft excess (whose origin is unclear) with 
respect to a Galactic absorbing column, more so in the pn; but that 
when the absorbing column is allowed to be free, the three detectors 
agree reasonably well with regards to the other cluster properties. 
When all three detectors are analysed simultaneously, the fits adopt 
a compromise between the pn and MOS values (there are roughly 
equal numbers of counts in the pn and two MOS combined). 

Model C has a fixed A'h, but freely varying redshift z. The 







A 


B 


C 




A'h 


2.49 


1 47+°-3 


2.49 


Ml 


kT 
Z 


o cc+0.05 
—0 05 
25+°'-°^ 

—0 02 


3.69+«-0^ 

o.26iu;»2 


T 55+0.05 
■ —0 05 

0-26+°:°! 




z 


0.0822 


0.0822 


n riTS+O-Ooi 






349.6/343 


340.5/342 


342.1/342 




A'h 


2.49 


1 91+°-3 


2.49 


M2 


kT 
Z 


, 43+0.05 
—0 05 
22^°'-°^ 

—0.02 


3 50+"-''^ 

—0 07 

23+°'-°^ 

—0.02 


3 42+0-OS 

■ —0 05 

22+°'-"^ 

■ —0 02 




z 


0.0822 


0.0822 


n 081+0-002 
'J-Uol -0.002 






297.3/348 


294.2/347 


296.8/347 




Nh 


2.49 


n 01+0-2 
"■'J^-o.oi 


2.49 


pn 


kT 
Z 


, ^^+0.030 
—0 003 
0.20+^-0; 

— U.Ul 


0.25+":"! 

— U.UJ 


n O9+0.02 
— 03 

o.22+«:"l 

—U.Ul 




z 


0.0822 


0.0822 


n 077+0-0004 

U.\J 1 1 _Q QQJ 






1321.2/1003 


1118.8/1002 


1293.6/1002 




Nh 


2.49 


0.59^»:1 


2.49 


All 


kT 
Z 


O TC+O.Ol 

-'■-'■^-0.005 
n 919+0.003 
"•^^^-0.006 


3.61^:2^ 

n 95+0-01 
"■^-"-o.oi 


, ,0+0.02 
-'--'-'-0.02 
n 93+0-01 
"■^^-0.01 




z 


0.0822 


0.0822 


079+0-001 




x' 


2000.6/1698 


1815.8/1697 


1968.0/1697 



Ah is Galactic column in units of 10^" cm"^- temperature in keV; Z 
metallicity relative to solar. Errors are 1 a. Quantities without errors were 
fixed. Models A-C are phabsx mekal, with: Galactic Ah; free A'h; free z. 
When fitting all three detectors simultaneously, the normalizations were 
untied (though the effects of this were slight in this instance). 

three detectors adopt disparate redshifts, with MOS 1 and the pn dis- 
agreeing (at the IcT level) with the nominal optical redshift, 0.0822. 
Analysis of the RGS data, however, with its higher spectral resolu- 
tion (see Section l?!2l . results in a redshift consistent with the opti- 
cal value. Comparing the results of models A and C (which differ 
only by whether redshift is free or not), we see that the only effects 
of freeing z are to reduce somewhat - the results for other pa- 
rameters such as T and Z are unchanged. Accordingly, we have left 
z fixed at the optical value in all subsequent EPIC fits. 

In Table |2| are shown the results of fitting the entire central 
4arcmin. Now that we include the cooling region, a wider variety 
of spectral models become relevant. Models A and B are the same 
as in the i-^' case. The fits in the 0-4arcmin region adopt lower 
temperatures and higher abundances, so that we may expect to see 
temperature and abundance profiles that decrease and increase re- 
spectively with decreasing radius (see Section|4}. 

The resu l ts of m odel B can be compared directly with those 
of lAllenetalJj200lh . who found, using A5CA: T = 3.38+o j],, Z = 
0.35+[5 []^, and A'h = 0.99+J] J58 (90 per cent confidence limits). 

In models D and E we add a second mekal component, with 
an independent temperature and normalization. For all three detec- 
tors, the fits are significantly improved. For example, for MOSl 
with a free A'h, the change in fit statistic is A;^^ = —17.9 for the 
removal of two degrees of freedom. According to the F-test, the 
significance of such an improvement is (1 — 3.4 x 10^^). Adhering 
to strict mathematical formalism, the F-test is not really valid in 
cases such as this (testing for the presence of a second component), 
where one of the hypotheses lies on the boundary of the parame- 
ter space l^rotassov et al. 2002). In practice, howe ver, Monte Carlo 
simul ations show that it gives acceptable answers jjohnstone et alj 
l2002t) . 
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Table 2. Global (0-4arcmin) spectral properties. 



A B D E F G 



Nn 2.49 1.49+5J-^ 2.49 2.52+jJ-^ 2.49 2.11+S5i 

kT 3 ^5+0.03 o 47+0.04 , 4f,+0.04 o 4<r+0.b4 , Sfi+O-O^ ^ "ifi+O BS 

HI j.jj_Qg3 •'■^'-0.04 -'■™-0.03 -'■^•'-0.03 -'-^"-0.05 •'■•^"-0.05 

Ml K^IK - - O.OlStlZ 0.015^«SJ«^ 

M - - - - 93li| 73^20 

7 n ^';+o.02 n ^7+0.02 n ^q+O-OI n ^q+0-02 r, ,7+0.02 r, 07+0.02 

u.jj "-^'-0.02 "--^'-0.01 "■-^'^-0.02 "--"-0.02 "-"-0.02 

X- 474.9/396 450.0/395 432.0/394 432.1/393 437.4/395 435.4/394 



Nn 2.49 l-^S^Ui 2.49 2.13^JJ-^ 2.49 2.56+°] 

IT T. 99 +0.03 o on+0-03 ^ 09+0.03 n o<;+0.b5 o ^ An+OOS 

(CJ j.z^ QQj j.ju_QQ3 j.j^ ,j(j3 j.jJ qq^ ■^■^"-0.04 -"-^"-0.04 

*:r2 - - 0.75+Ui, 0.86+02 

M2 K^IK - - 0.016^«:°«3 0.015+«:t 

M - - - - 92t\l 96+1 

7 n ^n+O-Oi n^i+O-Oi n ^9+0-02 n ^9+0-02 n^i+O-Oi mi+o^i 

A U.JU_Q(,[ U.JI_fl(„ U.JZ Q(,2 "■^''-0.02 "-^^-0.01 '^■^^-0.01 

404.3/400 391.9/399 375.5/398 374.0/397 372.0/399 371.9/398 



Wh 2.49 0-47iol 2.49 0.87+°-| 2.49 l-H^ol 

^T ^ 04+0-02 o 90+0.02 o 90+0.02 9 90+0.08 o /o+0.03 o 4^+0.03 

HI ^-"^-0.02 -'■^'-0.02 -'■^"-0.02 ^-^'-0.08 -'■^-'-0.04 ^-^-^-0.03 

m - - o.3i;|>:°i 6.78i« 

pn K./K - - 0.51+00^ 

M - - - - 159+*o 86+|] 

7 n 97S+0-008 n^ia+O-009 rn9';+0-0i 097+0.01 r,^ni+0-009 0^99+001 

;t;2 1781.8/1190 1482.2/1189 1537.3/1188 1307.4/1187 1504.2/1189 1417.2/1188 



Models: A, B phabs x mekal; D, E phabs x (mekal +mekal); F, G phabs x (mekal + mkcf low). The first, second model in each pair have Galactic, free 
Wh respectively. K^jK relative normahzation of second temperature component; M mass deposition rate ( yr" ' ); Z metalhcity (solar units). 

The relative normalization of the second temperature compo- 
nent is a few per cent of that of the bulk hot gas. If the hot and 
cold phases are in pressure equilibrium, then hqJc = nnTk- The 
normalization of the mekal model varies as = ii^V. Hence, the 
relative filling factor of the cold phase is given by 

Vc/Vh = (^c/^h) * {n^lncf = {KclKn) * {Tc/Tnf; 

which applied to the results in Table|2|for each detector gives values 
^ few X 10^'* in this case. 

Interestingly, when a second temperature component is avail- 
able, then the fitted Nh for the MOS detectors becomes consis- 
tent with the Galactic value. The pn results, however, for model 
E are not sensibly constrained. It is noticeable that the best-fitting 
low-temperature components for the two MOS detectors in mod- 
els D and E are somewhat discrepant. If (as we will argue later on) 
there really is a cooling flow operating down to low temperatures in 
this system, then there will be gas present at a continuous range of 
temperatures, and it is not obvious what the preferred single low- 
temperature value will be in a model with only two temperature 
components, nor how it will vary between different instruments. 
The upper temperatures for the two MOS instruments also seem to 
be systematically different, albeit at a lower level than for the low 
temperatures. In terms of the fit, the higher value for T2 in MOS2 
seems to be due to a slightly broader iron L hump at 1 keV in this 
detector. 

For completeness, we note that the pn results are stable against 
alternative forms of light-curve filtering and background correc- 
tion. For example, if we apply the same count-rate cut (0.45 cts^') 
as used in the blank-sky fields of Lumb ( 2002), and then employ 
a double background subtraction technique (as used for the MOS), 
then the pn fits are essentially unchanged. 

Using a cooling flow model (in which gas cools from the tem- 
perature of the bulk hot gas), instead of a single temperature, for 



the second component, also provides a significant improvement in 
fit over a one-component model. The data are not able to distin- 
guish between these two possibilities with any degree of accuracy. 
For example, for MOSl a cooling flow provides a very slightly 
worse fit than a single temperature, whereas for the MOS2 the sit- 
uation is reversed. The cooling flow fit does have only one extra 
free parameter though. When the absorbing column is allowed to 
be free, all three detectors independently fit a mass deposition rate 
~ 90 Mq yr^ ' . An analysis of the Reflection Grating Spectrome- 
ter (RGS) data (Section |5} produces very similar mass deposition 
results. 

In Fig.|3|we show the results of allowing the minimum tem- 
perature of the cooling flow component, Tnjin, to vary freely; for 
both fixed and free Nn (i.e. models F and G in Table|2|respectively). 
When Tjnin is allowed to vary freely, a value somewhat above 'zero' 
(i.e. the model minimum, 0.081 keV) is preferred; but a zero cut- 
off temperature is entirely consistent with the data, within the IcT 
limits. The associated mass-deposition rate, M, is little changed. 
Freeing Nn broadens the allowed confidence regions, but does not 
greatly affect the results. 



3.1 Comparison with Chandra 

We present a brief comparison of the results of XMM-Newton and 
Chandra observations of Abell 2597 through a re-analysis of the 
Chandra observation of McNamara et al. (2001). We use a con- 
servative flare screening of the Chandra data that results in 7ks 
of good time. The agreement between the Chandra and XMM- 
Newton pointings for the location of the X-ray emission peak is 
excellent. The largest complete circle centred on this point that 
can be extracted from the ACIS-S3 chip has a radius of 1 arcmin. 
This encompasses most of the cooling region (see Section 1421 . 
A background spectrum was extracted from the appropriate stan- 
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Figure 3. Confidence contours in the M-r,ni„ plane for the MOSl + MOS2 
coohng flow fits (models F and G in Tablelzl but with T^i^ free to vary). Up- 
per panel fixed Nh, lower panel free Nh. The cross shows the best-fit value, 
and the contours are plotted for Ax^ = 2.3, 4.61, 9.21; coiresponding to the 
68, 90, and 99 per cent confidence limits respectively on two free parame- 
ters. The left-hand boundary is due to the lowest usable model temperature, 
0.081 keV. 




0,5 



Figure 4. Confidence contours in the M-T^in plane for the central 1 arcmin 
in radius (the cooling region) from the Chandra ACIS-S3 data, with a 
Galactic Nu. Plot details are as per Figl3l 



dard blank-sky field. The spectrum was fitted over the 0.6-7.0keV 
range. 

Fig.Ois the Chandra counterpart of the upper panel (i.e. for a 
fixed Nh) of Fig.|3lfor the MOS. We were unable to obtain sensible 
constraints when the absorbing column was allowed to vary freely. 
The preferred low-temperature cut-off for the cooling flow compo- 
nent, though not tightly constrained, is as low as possible. The as- 
sociated mass deposition rate is --^ 70 Mp, yr^', entirely consistent 
with the XMM-Newton EPIC results for the central 4 arcmin. These 
results are also in agreement with those from the RGS, as shown in 
Fig. 1171 Fig|4|may be contrasted with the more common Chandra 
result in suc h cases, where a mu ch higher T^[^ is preferred (e.g. in 
Abell 3581: IJohnstone et all2005L fig. 10). 



4 RADIAL PROPERTIES 

Given the symmetric nature of the emission visible in Fig. |2| we 
have performed analyses using circular annuli centred on the clus- 
ter emission peak. 

4.1 Surface brightness profile 

Fig- El shows the surface brightness profile of the central 4 arcmin 
in radius of the cluster in the 0.5-2. OkeV band, from the sum of 
the MOSl and MOS2 data. The profile was extracted from 1 arcsec 
pixel images and vignetting-corrected exposure maps, using circu- 
lar annuli centred on the emission peak. A background was sub- 
tracted using the appropriate scaled blank-field data sets. The width 
of each annulus was such that the significance of the background- 
subtracted counts it contained was at least 5(T. Point sources were 
excluded. 

Also shown is the result of a j8 profile fit to the surface 
brightness, Z(r) = Zo[l + {r/r^f]-^P+^^-^- with j8 = 0.585+j];|][|^, 
= 23. Olojjj arcsec. 

The surface brightness profile is affected by the XMM-Newton 
PSF, which can itself be well described by a /3 model ( GhizzardJ 
2001 , 2002), with a core radius of around 5 arcsec. When an intrin- 
sic brightness profile a{7) is modified by an instrumental PSF b{r), 
the observed profile is the convolution of the two, 

c{r) = I (fua{u)b{r-u) =2k dk kJo{kr)a{k)b{k), (1) 

where the second form is for the case of radial symmetry and makes 
use of the convolution theorem. d( k) in this expressio n is the ra- 
dially symmetric Fourier transform IB irkinshawll 1 994^ , or Hankel 
transform of order zero, 

f{k)= r drMkr)rf{r), (2) 
■h) 

with ^0 a Bessel function of the first kind. Evaluating the effect of 
instrument PSF on a surface brightness profile therefore reduces to 
performing three Hankel transforms. The general /3 function has an 
analytic Hankel transform ( Gradshteyn & Rvzhik 2000, equation 
6.565.4), but the specific case of /3 = 2/3, which is a good approxi- 
mation for both the XMM-Newton PSF and many cluster brightness 
profiles, is especially simple: 

-ak 

f[r) = (a^ + r\^l^ ^ f{k) = 0) 

a 

As a result, if 

S{r) = So[l + (r/rs)2]-3/2 and P{r) = Po[i + {r/r^fr^^^ (4) 



describe the source and PSF profiles respectively, with the latter 
normalized so that = {Inrp ) ^ ' , then it is straightforward to show 
that the convolution c(r) is given by 




This is nothing more than another jS profile, with a reduced central 
amplitude 5q and increased core radius r'^ : 



S^, = 5()(l + ^) ~5o(l-2^) and ri = r, + rp. (6) 

The simplest PSF correction we can make to the surface 
brightness profile is therefore to subtract the radius of the PSF, 
5 arcsec, from the core fitted to the convolved profile, 23 arcsec. 
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Figure 5. Surface brightness profile of the central 240 arcsec in radius, 0.5- 
2.0keV band, from summed MOSl and M0S2 profiles. The dotted line 
shows a ^ profile fit to the data - see text for details. The profile has been 
corrected for background and vignetting. 



resulting in an estimate of 18 arcsec for the intrinsic surface bright- 
ness core radius. 



4.2 Spectral profiles 

In order to obtain reliable spectral fits, annular radii were chosen 
so as to encompass slightly more than 10000 (net, i.e. background- 
subtracted) counts in each MOS, and thus over 20000 counts in 
the pn. It was found that annuli with fewer counts (e.g. 7500 per 
MOS) tended to give poor (unphysical) results in the deprojection 
and PSF-correction analyses. In such cases, the lower signal-to- 
noise allows physically unrealistic solutions (in which fitted quan- 
tities tend to oscillate about the smoother profiles obtained from 
annuli with more counts) to become mathematically valid. Both 
correction algorithms (see e.g. jlohnstone et al. 2005 for testing of 
the deprojection process) involve redistribution of counts between 
annuli, potentially (depending on the details of the brightness pro- 
file) leaving small net numbers of counts in any one annulus. In 
practice it is found to be beneficial not to let the uncorrected count 
numbers fall too low. The outermost annulus was truncated at a ra- 
dius of 5 arcmin (due to the increasing relative background) and so 
contains somewhat fewer counts. By the time this last annulus is 
reached, the background is contributing about 30 per cent of the 
total counts. 



4.2.1 Projected profiles 

Fig. |5| shows the results of fitting a single-temperature MEKAL 
^Meweetal.lll995^ model to each spectral annulus, with an uni- 
form (but free to vary) absorbing column applied to all annuli. The 
central temperature is a little over 2.5 keV (a function of the size 
of the central annulus). Over the central lOOkpc or so the tempera- 
ture increases, reaching a maximum of ~ 3.6keV. Outside the cen- 
tral cooling region, there is perhaps evidence for a gradual decline 
in temperature beyond ~ 150kpc or so. The metallicity exhibits a 
fairly smooth decline with increasing radius in the central cooling 
region, dropping from around 0.5 Z0 at the centre to 0.2 Z,;, at the 
outside. 

Results are shown for the simultaneous fit of both MOS instru- 
ments, and also for the pn. The two instruments agree reasonably 
well in terms of T and Z, but have very different preferences for 
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Figure 6. Projected temperature (upper panel) and metallicity (lower panel) 
profiles. The open triangles are simultaneous MOSl + MOS2 fits, the filled 
squares are the pn. The model was phabs x mekal with a common, free 
absorbing column shared by all annuli. The Nw values were (in units of 
102°cm-2): 1.52+j5-558 (MOS); 0.19+5J o6 (pn). WithA^H free to vary between 
annuli, the T and Z profiles are not much changed. The MOS A'h shows no 
significant variation with radius; whereas that of the pn is consistent with 
the MOS in the central two bins, but much lower outside this region. 



In Fig. we show the radial dependence of the cumula- 
tive mass deposition rate, M(< r), within the cooling flow radius 
(130kpc; see Section l4.2.2t . There is a discrepancy between the 
MOS and pn detectors, with the former preferring smaller mass 
deposition rates. Owing to the way this figure was constructed (fit- 
ting each annulus independently with a single temperature compo- 
nent and a cooling flow component, then summing the M values so 
obtained from the centre outwards), the discrepancy between the 
two instruments buil ds with increasing ra dius. Also shown is the 
M value inferred by lOeeerle et alJ i200ll) from their FUSE data; 
namely 40 M© yr^' within a radius of 40 kpc. The figure illustrates 
that there is reasonable agreement between the X-ray and UV mass 
deposition rates within a radius of 40 kpc. 



4.2.2 Deprojected profiles 

The results of Fig. |6| are subject to projection effects, in which 
the spectral properties at any point in the cluster are the emission- 
weighted superposition of radiation originating at all points along 
the line of sight through the cluster. In Fig.|8|we show the results of 
a deprojection of the MOS and pn spectra using the XSPEC pro j ct 
model. Under the assumption of intrinsic spherical (more generally, 
ellipsoidal) shells of emission, this model calculates the geomet- 
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Figure 7. Cumulative projected mass deposition rate, M(< r), for the cool- 
ing flow region. The open triangles are simultaneous MOS 1 + MOS2 fits, 
the filled squares are the pn. The model was phabs x (mekal + mkcf low) 
w ith a fixed Galactic ab sorption. The star shows the FUSE O VI data point 
of lOegerle et alJ 1200 ih . with an error bar appropriate for the accuracy of 
the FJ/Sii flux measurement, (1.3±0.35) x 10~'^ ergcm"- s^'. 



ric weighting factors according to which emission is redistributed 
amongst the projected annuli. 

As is to be expected, the temperature in the outer regions 
where the profile is fairly flat is seen to be relatively little influenced 
by projection. Deprojection recovers a lower central temperature 
than before, since in the projected fits the spectrum of the central 
annulus is contaminated by overlying hotter emission. There is also 
some evidence for a higher central abundance in the deprojected 
metallicity profile. 

The deprojected pn temperature profile shows some signs of 
instability in the centre, in which neighbouring bins 'oscillate' with 
high and low temperatures. The interpolated values, though, agree 
well with the MOS results. Recently, Johnstone et al. ( 2005) have 
shown that the pro j ct approach works well at reproducing various 
simple synthetic cluster temperature and density structures. Possi- 
bly the slight instability of the pn fit reflects a more complex under- 
lying temperature distribution (e.g. more than one phase); but more 
likely it is simply a quirk of geometry that renders certain phys- 
ically unrealistic solutions mathematically plausible. Freezing the 
abundances at the projected value makes little difference, but the ef- 
fect can be removed by using larger annuli (although the pn annuli 
shown already contain similar numbers of counts to the combined 
MOS annuli), but we retain the same annuli here for both instru- 
ments for illustration. 

In Fig.|9|we show various quantities derived from the depro- 
jected spectral fits. The electron density «e is obtained from the 

1/2 

mekal normalization, and scales as 

7-1/2, 



jQ . The cooling time (scal- 
ing as Hq ''^) is an approximate isobaric one, calculated as the 
time taken for the gas to radiate its enthalpy using the instantaneous 
cooling rate at any temperature, 

kT 



^cool ' 



H 



Y 



«e«HA(r) y- 1 /XXH«eA(r) 



(7) 



where: y = 5/3 is the adiabatic index; /i ~ 0.61 (for a fully-ionized 
0.3 Zq plasma) is the molecular weight; Xn ~ 0.71 is the hydro- 
gen mass fraction; and A(r) is the cooling function. The 'entropy' 



is the standard astronomer's entropy, 5 = Taj, 



-2/3 



and scales as 
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Figure 8. Deprojected temperature (upper panel) and metallicity (lower 
panel) profiles. The open triangles are simultaneous MOS 1 -I- MOS2 fits, the 
filled squares are the pn (as per Fig.|6|. The model was proj ct x phabs x 
mekal with a common, free absorbing column shared by all annuli. The 



Nn values were (in units of 10^" cm"^): 
little changed from the projected results. 



, ,ro-l0.08 

^••'-'-0.08 



(MOS); 0.23 



4-0.06 



(pn); 



(in perfect agreement with the ROSAT result o f ISarazin et aljl995l) . 
and less than 5Gyr inside 85 kpc. 

Also shown in Fig.|9|are the results of simple model fits to the 
density and entropy profiles. The density profile is represented by 
a j3 model, n(r) = no[l + (r/rc)^]-'-^'^; with /3 = 0.57 



rc = 49.0i'iy n kpc. These results for the den- 



-1-0.01 
-0.01 



no ' 



03^4+0.0006 -3 
"•^•'■'^-0.0006 ^™ ' 
sity profile are in reasonable agreement with those obtained from 

fitting the surface brightness profile (Section f4.U . The absence of a 
significant excess in the deprojected density in the last annulus in- 
dicates that we reach a sufficiently large radius that projection from 
more distant material is not significant (compa re with the Chan- 
dra deprojection results of ljohnstone et alj2005h . The entropy pro- 
file is parametrized as a constant plus power-law, S{r) = So + kr^ ; 
with 5o 



15.1 + [-^keVcm2, 



yt = 0.28+!!l]^keVcm2kpc 



-1.19 



The cooling time is less than lOGyr inside a radius of 130 kpc 



1.19_Q Q^. The fit with the 5o term is substantially better than that 
without, but bear in mind that the PSF has not been corrected for 
here. A logarithmic slope for the entropy profile of ~ 1.1 is pro- 
duced i n simulations involving gravitational collapse and shock 
heating iTozzi & NormarJ2'00U) . 

The virial radius becomes a quantity of interest at this point. 
It is known that the calibrations of the scaling relations obtained 
from simulations have normalizations quite different to those ob- 
tained from observation. Allen et al. 1 200 IJ studied the cluster 
virial scaling relations using Chandra, but do not quote an explicit 
radius-temperature relation. Using their data for an SCDM cos- 
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Figure 9. Derived quantities from the deprojected profiles of Fig. |8] up- 
per panel electron density; middle panel approximate isobaric cooling time; 
lower panel entropy. The dotted lines show the fits of some simple models 
to the data: a beta profile for the density; and a constant plus power-law for 
the entropy. See text for details. 



mology, we calibrate the r-T relation as ?'2goo = 0.31Mpc /j^q (1 + 
^) -3/2 (jjjQQ / keV)^'^ . The conversion from an overdensity of 2500 
to one of 200 (appropriate for the virial radius) depends upon the 
form assumed for the density profile. Taking an NFW profile with 
a typical cluster concentration c = 5, we find that r2oo ~ 2.75^2500. 
Using a temperature of 3.5 keV, we therefore estimate a virial ra- 
dius of 1.4Mpc for Abell 2597. The radius 0.1i?v therefore lies 
within the sixth annulus from the centre (a bin whose mid-point is 
at 150 kpc). The measured entropy in this bin is 5o ir-^/ = (137 ± 



6)/i,n'^ keV c m . T his is a relatively low value (see e.g. fig. 4 of 
IPoimian et alJl20()3l) . but recall PSF effects are still present in this 



result. 



4.2.3 PSF-corrected profiles 

The broadening effect of the XMM-Newton PSF acts to redistribute 
counts between spectral annuli in a manner conceptually similar to 
that of projection. Consequently, it can be corrected for in the same 
way, namely by using an XSPEC mixing model, ximnpsf . 

The most recent version available at the time of writing was 
unable to fit observations from different instruments simultane- 
ously, so we present here only the results for MOSl, in order to il- 
lustrate the effects of PSF on the spectral fits. We have been unable 
to achieve physically realistic solutions when using the xmmpsf and 
pro j ct models in combination, to correct for both projection and 
PSF at the same time. The exception is when annuli that are phys- 




Region 

Figure 10. xmmpsf mixing factors calculated at l.OkeV. Each curve shows, 
for one of the annular regions, the fraction of the flux in that region origi- 
nating/i-om each of the other regions. Only contributions greater than 1 per 
cent are shown. The peak flux is always from the region itself. The actual 
mixing factors used internally by the model relate to the fraction of the flux 
from each region going to other regions. The factors shown here are calcu- 
lated from these using the model normalizations fitted to each annulus. By 
construction, the points with the same x-coordinate sum to 1.0. 



ically very large are employed; but in such cases the PSF redistri- 
bution becomes less significant anyway, and we also obtain very 
limited resolution of the central cooling region, which is the main 
area of interest here. 

For a given set of spectral annuli, and a given brightness pro- 
file, the xiranpsf model calculates the mixing factors through which 
flux from any given annulus is redistributed to all the others. The 
factors are calculated at ten energies in the range 0.5-6.0keV, al- 
though the energy dependence of the XMM-Newton PSF is rela- 
tively small. The input brightness profile can be specified using ei- 
ther an image, or a model parametrization. We find that there is no 
significant difference between using a MOSl (say) image or the 
PSF-corrected /3 profile from Section im to provide the brightness 
information. By way of example, we show the mixing factors for 
our set of annuli at l.OkeV in Fig. llOl At the ten per cent level, each 
annulus is only affected by its immediate neighbours. Nevertheless, 
the mixing is potentially quite significant, since for many regions 
only 50 per cent of their flux originates from that same region. 

In Fig.ET] we show the T and Z profiles obtained from a PSF- 
corrected fit of the MOS 1 data. Compared to the uncorrected pro- 
files (Fig.|6j, a steeper central temperature gradient is recovered, as 
well as a larger central abundance peak. 

Fig. ll2l displavs the density, cooling time, and entropy derived 
from the PSF-corrected profiles. Compared to the results of Fig.|9| a 
higher central density, and a lower central cooling time and entropy 
are obtained. Also shown is a jS model fit for the density, with: j8 = 
0.579+0 004, no = 0.068lo;oo} cm-^, = 31.6+o;^kpc. The PSF 
corrected core radius of 15arcsec is comparable to the estimate of 
ISarcsec made in Sectionj^Jusing the surface brightness profile. 

The entropy profile is well-fit - reduced = 7.49/(9 — 3) 
- by a constant plus power-law model of the same form 
used in Section 14.2.21 with: Sq = 7.9+} jkeVcm^, k = 
0;08keVcm2kpc-'i5 



0.40^ 



, 77 = 1.15 



hO.OS 
-0.04- 



4.3 Comparison with Chandra temperature profile 

As per Section ism we have compared the XMM-Newton and Chan- 
dra results. Chandra spectra were extracted from the same annuli 
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Radius / kpc 

Figure 11. PSF-corrected MOSl temperature (upper panel) and metal- 
licity (lower panel) profiles. The model was xmmpsf x phabs x mekal 
with a common, free absorbing column shai'ed by all annuli: Wh = 
1.40^0;] 1020 cra-2. 



as used in the XMM-Newton analysis. Restricting ourselves to those 
annuli that lie entirely on the ACIS-S3 chip, we obtain five spectral 
regions from the Chandra data. Background spectra were produced 
as previously described (though background should not be signifi- 
cant in this region). 

In Fig. ll3l we compare the temperature profiles obtained from 
MOSl and Chandra, using a simple phabs x mekal model with 
a fixed Galactic absorption. The agreement between the MOSl 
and ACIS temperatures is good everywhere except in the central 
bin, where the ACIS temperature is significantly lower than that of 
the MOS. Correcting for the XMM-Newton PSF using the xmmpsf 
model in the previously described manner lowers the central MOS 1 
temperature so that it agrees very well with that of Chandra. The 
agreement of the other MOS 1 points with those of ACIS is wors- 
ened though, as the PSF correction raises all but the central MOS 1 
temperatures. 

The XMM- NewtonI Chandra cr oss-calibration has recently 
been examined bv lKirschetai]<2004 . 
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Figure 12. Derived quantities from the MOSl PSF-corrected profiles of 
Fig.[n] Details as per Fig. 191 
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Figure 13. Comparison of MOSl and Chandra ACIS-S3 temperature pro- 
files. Open symbols MOSl: circles non-PSF corrected; stars PSF corrected. 
Filled squares Chandra. The model in each case was phabs x mekal with 
a fixed Galactic Nw . The PSF connection was achieved with the addition of 
an xmmpsf component (see text for details). 



5 REFLECTION GRATING SPECTROMETER 

Around half of the light in each telescope feeding the MOS de- 
tectors is diverted to a Reflection Grating Spectrometer (RGS; 

III I 
den Herder et al. 2001). Each RGS contains 9 CCDs (though one 

has failed in each RGS) along the dispersion direction, and is sen- 
sitive to soft X-rays in the range ~ 6-38 A. The field of view in the 
cross-dispersion direction is 5 arcmin. 



5.1 Data reduction 

The two RGS were operated in the standard Spectroscopy mode. 
Data were processed from the ODFs using the SAS meta-task 
rgsproc. The position of the cluster emission peak as measured 
from the MOS image was used to specify the source extraction 
point. The generated RGSl and RGS2 events files were filtered 
for periods of high background with the same approach as used i n 
the construction of the RGS background files iTamura et all2003l) . 
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Figure 14. Profile of RGSl counts in 10 arcsec bins in the cross-dispersion 
direction, for the events lying in the first order PI-j3 (dispersion angle) se- 
lection region. 



Indeed, since the background files do not have the necessary in- 
formation for further filtering, we have little choice but to proceed 
in this manner (though we might prefer to use a stricter filtering in 
this case). That is, we form a light-curve from the (FLAG == 8,16) 
events at high cross-dispersion angles (|XDSP_CDRR| > 1 . 5e-4) on 
CCD9 (that which covers the highest energy band and is most af- 
fected by background), using 100 second time bins. Only periods 
with count-rates less than 0.15 counts per second for each RGS 
were retained, leaving ~ 74 ks of time per RGS. Observe that al- 
though this is significantly larger than the remaining EPIC expo- 
sure, both EPIC and RGS have (essentially) been filtered to the 
same level as the appropriate ESA-supplied background templates 
(soft protons channelled through the mirrors are responsible for 
most of the XMM-Newton background flaring, and relatively few 
are scattered by the gratings into the RGS detectors). 

Since the field of view in the cross-dispersion direction is only 
5 arcmin wide, it is essentially filled by cluster emission from Abell 
2597. We theref ore make use of the standard background template 
files provided bv lTamur~e t al. ( 2003). We note that these files were 
made at an operating temperature of -80°C, but that in late 2002 
both RGS were cooled to -1 10°C. 

In Fig.fT4l we show the 'one-dimensional surface brightness' 
profile obtained from RGSl by summing counts in 10 arcsec bins 
in the cross-dispersion direction. The centre of Abell 2597 is un- 
fortunately offset from the mid-point of the CCD array by about 
1 arcmin. There are insufficient counts to form profiles in narrow 
energy bands, in order to examine the possible differences between 
the profiles of differen t spectral lines, e.g. due to resonant scattering 
jSakelUou et all2002h . 

RGS spectra are extracted by making selections both in 
the spatial (dispersion/cross-dispersion) and energy (dispersion/PI) 
planes, by making use of the intrinsic energy discrimination of the 
detector CCDs. Due to the extended nature of the source, we have 
used slightly broader selection criteria than the defaults. Our spatial 
selection took 97 per cent of the PSF in the cross-dispersion direc- 
tion, which is equivalent to the central 2 arcmin or so of the cluster, 
i.e. the central cooling region. In the energy plane we took 94 per 
cent of the pulse-height distribution. We restricted our analysis to 
the first order spectra. 

Response files were generated using the SAS task rgsrmf gen 
with standard parameters (4000 energy bins between 0.3 and 
2.8keV). Spectra were grouped to a minimum of 20 counts per bin. 
RGS source and background spectral channels can have different 



qualities, and since XSPEC ignores quality information in the back- 
ground spectrum, the source PHA files were modified by adding in 
any extra bad channels from the background spectra. We fit the two 
RGS instruments simultaneously over the 5-38 A range, thereby 
compensating for the missing iron L complex and O VIII line ab- 
sent (due to failed CCDs) from RGS 1 and RGS2 respectively. 

In an attempt to deal with the spectral blurring that results 
from the extended nature of the source, we have employed the 
XSPEC rgsxsrc model (due to A. Rasmussen). This convolves the 
spectral model with an angular structure function computed from 
an image (we used the MOS 1 image in the RGS energy range of 
0.3-2.5 keV) of the target region. Correcting for broadening in this 
way is only an approximation to the complex instrument response 
to a spatially extended input. The rgsxsrc model assumes that the 
spatial structure of the source is independent of energy. To test the 
possible consequences of this assumption, we examined the effects 
of using input images in different energy ranges (e.g. below IkeV, 
or dividing the RGS energy range in two and using either the up- 
per or lower half). The results quoted in subsequent sections were 
found to have very little dependence on the energy range used for 
the rgsxsrc image. For example, the changes in fitted fluxes for 
the FeXVII lines (see Section l531 when using a different energy 
range for the rgsxsrc image are typically of order a few per cent, 
in a random direction (i.e. not systematically higher or lower for an 
image in any given energy range). The changes in line significances 
were negligible. 

According to the XMM-Newton Users' Handbook^ , the the- 
oretical wavelength resolution of the RGS spectral order n for an 
extended source of angular size 6 (arcmin) is degraded accord- 
ing to the formula Al/A ^ 0.138 B /n. Thus, for a source region 
with extent ~ 2 arcmin we may expect a wavelength resolution 
~ 0.25 A. This is borne out by examination of the model width of 
an rgsxsrc-blurred line of zero intrinsic width. 

5.2 Broad-band spectral fits 

In Fig. Us] we show the best-fitting single temperature vmekal 
model and its residuals as applied to the data. The abundances of 
iron and oxygen were allowed to vary independently (both these el- 
ements exhibit strong line features in the RGS energy range, so we 
may reasonably expect to obtain reliable constraints on their indi- 
vidual abundances), whilst those of all the other elements were tied 
together. The redshift and absorbing hydrogen column were also 
allowed to vary. Good fits (even with a single temperature model) 
could not be obtained when the absorption was fixed at Galactic. 
The fit parameters are listed in Table|3| The preferred absorption is 
substantially higher than Galactic (and that favoured by the EPIC 
instruments), for reasons that are not clear. The excess of above 
Galactic is consistent with the intrinsic A'h (Section^ detected by 
[b'Dea et al. 1 1994); but that was only seen against the small cen- 
tral radio source, so this is probably not significant. Allowing 
to vary results in temperatures and abundances that are consistent 
with those of the EPIC fits to the similar spatial region. 

5.2.1 Redshift issues 

Two simila r, but distinct, redshifts are in common use for 
Abell 2597.lKowalski eta'i] il983h . as used in the compilation of 
IStruble & Roo j ^99^) ~lneasured the optical redshifts of three 



|http : //xmm. vllspa. esa. es/external/xmm_user_support/docuinentatlon/uhb/ 
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Figure 15. RGSl (dark) and RGS2 (light) spectra, in units of ctA"' s"'. 
Also shown is the best-fitting single temperature model, together with its 
residuals (lower panel). 



Table 3. Results of RGS fits to the central ^ 2 arcmin region, 5-38 A. 
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Ny{ is Galactic column in units of 10^° cm^^; kT temperature in keV; Z 
metalHcity relative to solar for O, Fe, and all other elements (x); and z 
redshift. Errors are \c. All models are rgsxsrc x phabs x (M), where M 
is: A vmekal; B vmekal + vmekal; C -vmekal + vmcf low (with 
■vmcf low: T,^[ix = vmekal:?", vmcf low:rn]i„ = 0.081); D as C but with 
Tmm free- Where two phases were used, the metallicities were tied. 

(non-cD) cluster galaxies as z = 0.0874,0.0832,0.0851, hence z = 
0.0852 ± 0.002 (relative to Local Group standard of rest). Convert- 
ing this result to heliocentric using the NASA/IPAC Extragalactic 
Database (NED) velocity converter'* subtracts ~ 140 km s~' . givin g 
Zhelio = 0.0847. In contrast. lNoonanHl98lh . using lSchmidiHl965h . 
has z = 0.0821 (h eliocentric), based on the central radio galaxy. 
lOwen et alj jl995h find a heliocentric z = 0.0822 ± 0.0002, using 
optical emission lines of the central radio galaxy PKS llill-ll. 
The best-fit redshift found by Voit & Donahue ( 1997) for the opti- 
cal emission line nebula in Abell 2597 is z = 0.0821 ± 0.0002. In 
light of all this, we choose to adopt z = 0.0822 for the heliocentric 
redshift of the cluster. 

The best-fitting RGS broad-band redshift is consistent with 
the adopted optical redshift. The most prominent single line in the 
spectrum is that of O VIII Ka at around 20.5 A. The rest energy of 
this feature is 653.6eV. Fitting the spectrum in the 19-22A range 
with an rgsxsrc-blurred power-law plus a (zero intrinsic width) 
Gaussian, we find the best-fit value for the energy of this line to be 

"* |http : //nedwwH. ipac ■ caltech. edu/f orms/vel_correction.htmr| 
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Figure 16. Contributions of the two temperature components (solid lines) of 
model B in Tablel3lto the total (dotted line) model. The minority component 
is the low-temperature one. 

E = 603.8+[];^eV, implying a redshift of z = 0.0825+[J;i]J]jJ^, con- 
sistent with the broad-band RGS fit, and the optical value. 

We have also fit just the iron L complex, by restricting atten- 
tion to the 9.2-17A range (chosen to exclude the Mg line around 
9 A, and to include as much line-free continuum around the Fe 
L complex as possible). The plasma properties are not well con- 
strained in such a narrow wavelength range, with the exception 
of redshift, which is our sole interest here. According, we froze 
A'h at the broad-band value, and fit the data with a single mekal 
model. The best-fitti ng redshift for the Fe L complex was z = 
0.0831+0 001^. As per lStill & Mushotzkvli20o3) . the best-fit O VIII 
redshift is formally less than that of the Fe L complex. Within the 
1 CT limits, however, the redshifts obtained from the various RGS fits 
(broad-band, Fe L complex, O VIII) are consistent, and agree with 
our adopted optical reds hift. The hypothesis of s eparate velocity 
components in the core iStill & Mushotzkvll2002h could therefore 
be supported by these data, but does not seem to be required. 

5.2.2 Additional fit components 

Also shown in Table |3| are the results of adding a second temper- 
ature component to the fit, with an independent temperature and 
normalization. The removal of two degrees of freedom results in a 
change in fit of Ax^ = —28.5. According to an F-test, this improve- 
ment in fit is significant at the (1 - 3.1 x 10^^) level. In Fig. 1161 
we show the contributions made by each of the two temperature 
components to the total model. The bulk of the spectral fit is ob- 
viously dominated by the high-temperature component. The low- 
temperature component contributes several emission lines at wave- 
lengths around 15 A - these are examined in more detail below. 

In Table |3| we also show the results of fits using a cooling 
flow (vmcf low model) rather than a single temperature for the sec- 
ond component. These also provide significant improvements to the 
single-temperature fit, although the data are not of sufficient qual- 
ity to unambiguously state that a cooling-flow fit is preferred. The 
most interesting aspect of the fit is the obtained mass deposition 
rate: M ~ 100 Mg, yr"' . This is consistent with the EPIC fits to the 
central regions, Table|2| 

Allowing the low-temperature cut-off of the cooling flow com- 
ponent, r,^in, to be free (rather than fixed at the lowest value pos- 
sible in the XSPEC model, 0.08 IkeV), makes no improvement to 
the fit quality. Indeed, as we show in Fig. lI7l a very low temper- 
ature cut-off for the cooling flow is preferred. The slight upward 
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Figure 17. Confidence contours for the RGS cooHng flow fit (model D in 
Table l3l in the M-7^„in plane. Plot details are as per Fig|3] Compare also 
with the Chandra results of FigRI 

curve of the significance contours does indicate that higher mass 
deposition rates are permissible (though not preferred) with higher 
low-temperature cut-offs. 

5.3 Fe XVII line fitting 

We now return to the low-temperature emission lines in the 15 A re- 
gion of the spectrum previously alluded to. The Fe XVII lines in this 
region are particularly important because they are strong indicators 
of gas at temperatures ~ 0.3 keV, and consequently their presence 
is an important prediction of the cooling flow model. 

In order to examine these lines in more detail, we fit the 15- 
20 A region with a blurred power-law and four (zero intrinsic width) 
Gaussians. Adopting a fixed redshift of z = 0.0822, the line ener- 
gies were fixed at the positions appropriate for the 15.02, 16.78, 
17.08 A (rest wavelength) FeXVII fines, and the 16.00 A O VIII K/3 
line. We found the best-fit and IcT limits on the normalizations of 
each line. We obtained positive (but weak) evidence for all three 
important Fe XVII lines, at the l~2cT level (too low to qualify for a 
formal detection). 

The normalizations were converted to line fluxes (using a lu- 
minosity distance of 503 Mpc) by subtracting off the continuum 
component in the absence of a line. Results are shown in Table|4| 

We have also converted these fluxes to equivalent cool- 
ing flow mass deposition rates. This was achieved by using the 
mkcflow model in XSPEC to make a fake spectrum for a fidu- 
cial 1000 M0yr^' cooling flow, at a redshift etc. appropriate to 
Abell 2597. We fit the 15-20 A range of this model spectrum in the 
same way as applied to the actual data, in order to obtain the flux 
expected in each of the Fe XVII lines. The ratios of the Abell 2597 
line fluxes to those from the model were used to estimate the asso- 
ciated mass-flow rates. The results, as shown in Table|4| are poorly 
constrained, but consistent with both the broad-band RGS fits, and 
the EPIC fits to the central region. 

Table |4| also gives the F-test significance (the comments of 
Section |3] on the F-test also apply here) for each of the Fe lines. 
This is obtained from the improvement in fit that results when fit- 
ting a model with just the O VIII and each Fe line in turn present, 
compared to the best-fit with just the O line. The O line has been 
treated separately due to its different temperature dependence. The 
significance of each individual Fe line is admittedly low, varying 
from about 1-2(7. In order to try and achieve greater significance, 
we have tied the relative normalizations of the three Fe XVII lines 



Table 4. RGS Fe xvil line fits. 



A/A 


Power / lO'^'^ergs ' 


M/ M:.,yr-' 
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15.02 
16.78 
17.08 
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14+15 

19+12 
1^-14 
23+12 


34+" 
100+5] 
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Column 1: rest wavelength; column 2: power in line (Dl = 503 Mpc); 
column 3: mass-flow rate in line, inferred by comparison with a fiducial 
isobaric cooling flow model (see text for details); column 4: F-test 
significance (see text). 
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Figure 18. Ratio of data to a simple power-law model (i.e. no emission 
fines) in the 15-20 A region. The data have been re-binned, and the two 
RGS instruments have been combined (this requires a degree of interpo- 
lation since the wavelengths bins in the two instruments do not coincide 
exactly), for display purposes only. The sohd line shows the ratio appropri- 
ate for the Gaussian emission lines fitted in Table|3 the dotted line that for 
the cooling flow model C of Table l3l Note that the models were fitted to 
the un-binned data for the two separate instruments simultaneously (i.e. the 
models were not fitted to the data as displayed here). 



(all with similar temperature dependences) to that present in the 
fiducial isobaric cooling flow model, and just allowed the overall 
normalization to vary. The results are given in the 'All' row of Ta- 
ble|4| When combining the lines in this way, the significance is still 
low, but exceeds 2ct. 

In Fie. llSl we show the 15-20A region of the RGS spectrum, 
in terms of the ratio of the data to a model with no emission lines. 
Also shown are the appropriate curves for the model fitted for Ta- 
ble|4](power-law plus Gaussians) and for model C of Table |3| (iso- 
baric cooling flow). We note that the normalizations adopted by the 
Gaussians (where the line strengths may vary independently) are 
very similar to those of the cooling flow model. The exception is 
the 15.02 line, which is the Fe XVII line most susceptible to reso- 
nant scattering (e.g. Rugge & McKenzi3l 19851 see also the APED 
database^), for which a weaker fit is preferred. The quality of the 
line fits is clearly poor, but the values obtained for mass deposition 
rates agree with those from broad-band fits in both EPIC and the 
RGS. 



^ |http : //cxc . harvard ■ edu/atomdb/ 1 
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6 DISCUSSION 

First, we must caution tliat tiie overall quality of our observation 
is low, and that the formal statistical significance of our results is 
not large. In particular, the results on the FeXVII lines are really 
no more than suggestive. Nevertheless, the data do present a self- 
consistent picture between the various detectors, with an extremely 
interesting interpretation. 

Fitting each of the three EPIC detectors individually for a large 
central region covering that radial range where the cluster emission 
is high, then the addition of a cooling flow component provides 
a highly significant improvement in fit over a single temperature 
model (although we cannot statistically distinguish a cooling flow 
component, with a range of temperatures, from a single temperature 
second component). The best-fitting mass deposition rate is around 
9OM0yr"', with a la error of about 15 Moyr"' (see Table|2|for 
details). The cooling time is less than lOGyr within a radius of 
130kpc (Fig.|9}. 

Very similar results are obtained from the broad-band simul- 
taneous fit to the two RGS detectors (Table|3j. The main evidence 
in the data for low-temperature gas comes from the spectral region 
around 15 A, in particular the emission lines of FeXVII at 15, 17 A 
rest wavelength. We obtain weak (1-2(7; below formal detection 
levels) evidence for the three most important lines, and measure 
the power associated with each line (see Table |4}. By comparing 
the line fluxes with those produced from a classical isobaric cool- 
ing flow model, we can estimate the corresponding mass-flow rate 
for each line. Once again, we obtain values that are consistent with 
the 90 Mp; yr^' mark. Fixing the relative strengths of the three Fe 
lines to that predicted by the standard isobaric cooling flow model, 
and allowing the absolute normalization to vary, we obtain a fit sug- 
gesting the presence of Fe XVII at a significance of just over 2cT. 

We therefore establish, using three distinct methods, evidence 
for a cooling flow at a level --^ 90 yr^ ' . The preferred low- 
temperature cut-off for the flow, established from the RGS data, 
is essentially zero (0.081 keV, the minimum temperature available 
in the XSPEC mekal model). A plot of the M-T^jn confidence con- 
tours is shown in Fig. 1171 

A re-analysis of the Chandra data for Abell 2597 
iMcNamara et all200lll . as described in Section lTTI also supports 
the existence of a cooling flow at these levels, with a very small 
low-temperature cut-off - see Fig.|4| 

This ma s s depo sition rate is consistent with the results of 
lOegerle et alJ i200lh . using FUSE UV observations. These au- 
thors detected the O VI 1032 A resonance line (characteristic of gas 
at temperatures ~ 3 x 10^ K) in Abell 2597. Converting the de- 
tected flux to an equivalent mass deposition rate gives a UV mass- 
deposition rate ~ 40 Mp^ yr^' within the FUSE effective radius of 
40kpc. This is consistent with the X-ray results (Fig.0. 

In contr ast, Ovi was not detected in Abell 1795 

toegerle et al. 200 1|), nor in A bell 2029 or AbeU 3112 
tLecavelier des Etangs et alj |20o3)- Both from an X-ray and 
a UV standpoint, therefore, Abell 2597 appears to be an atypical 
object. The combination of the XMM-Newton X-ray and FUSE 
UV results for Abell 2597 suggests that it may harbour a classical 
cooling flow in which gas cools from ~ 4keV by at least two or- 
ders of magnitude in te mperature. Interestin gly, a recent Chandra 
analysis of Abell 2029 <Clarke et ai]l2004 shows that the X-ray 
data in this object are also consistent with a modest (by traditional 
standards) cooling flow extending down to very low temperatures, 
despite the lack of a UV detection. 

Among the many (e.g. iFabian et alj 1200 it IPeterson et alJ 



l200lt) explanations suggested for the 'cooling-flow problem', heat- 
ing from a central AGN, mediated by the rise of b uoyant plasma 
bubbles thr ough the ICM, is a pop ular candidate fe.g.l Vecchia et alJ 
2004; Rev nolds et ai]l2005t iRuszkowski et al. 2004). The Chan- 
dra observation of Abell 2597 iMcNamaraet al. 2001) showed it 
to contain X-ray surface brightness depressions, interpreted (as a 
result of the extension of spurs of old, low-frequency radio emis- 
sion) as 'ghost cavities', associated with a radio outburst ~ 10^ yr 
ago. 

Abell 2597 therefore appears to harbour a central AGN that 
generates buoyant bubbles, and yet it also appears to contain gas 
cooling to very low temperatures. The process of bubble genera- 
tion is necessarily an episodic one, and it therefore perhaps un- 
reasonable to assume that a steady, time-independent state will be 
maintained. Clusters may go through cycles of behaviour, in which 
cooling builds up to some level, then an AGN outburst initiates a 
heating phase in which cooling is suppressed (else massive cooling 
flows would be common phenomena). Abell 2597 would then be an 
object near the 'cooling catastrophe' point of the cycle. The relative 
frequency of such objects would depend on the relative timescales 
of the cooling and heating phases. The fact that such objects seem 
to be quite rare indicates that the phase of strong cooling is short- 
lived. 

Such a mechanism would require coupling between the state 
of the ICM and the activity of the central AGN, in order that a 
feedback loop could be maintained. 
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